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Abstract 

We consider the wide class of few-particle systems that have some analog of the thermo- 
dynamic laws. These systems are characterized by the distributions that are determined 
by the Hamiltonian and satisfy the Liouville equation. Few-particle systems of this class 
are described by a non-holonomic constraint: the power of non-potential forces is directly 



proportional to the velocity of the elementary phase volume change. The coefficient of 
this proportionality is determined by the Hamiltonian. In the general case, the examples 
of the few-particle systems of this class are the constant temperature systems, canonical- 
dissipative systems, and Fermi-Bose classical systems. 



> ' 1 Introduction 

in 

IT) 

■^j- . The main aim of statistical thermodynamics is to derive the thermodynamic properties of 

systems starting from a description of the motion of the particles. Statistical thermodynamics 

of few-particle systems have recently been employed to study a wide variety of problems in the 

field of molecular dynamics. In molecular dynamics calculations, few-particle systems can be 

exploited to generate statistical ensembles as the canonical, isothermal-isobaric and isokinetic 

ensembles [1, 2, 3, 4, 5, 6, 7, 8]. The aim of this work is the extension of the statistical 

thermodynamics to a wide class of few-particle systems. We can point out some few-particle 

H ' systems that have some analog of thermodynamics laws. 

c3 ; 

• The constant temperature systems with minimal Gaussian constraint are considered in 
Ref. [1, 2, 3, 4, 7]. These systems are the few-particle systems that are defined by the 
non-potential forces in the form F- n ^ = — ypi and the Gaussian non-holonomic constraint. 
This constraint can be represented as an addition term to the non-potential force. 

• The canonical distribution can be derived as a stationary solution of the Liouville equa- 
tion for a wide class of few-particle system [12]. This class is defined by a very simple 
condition for the non-potential forces: the power of the non-potential forces must be di- 
rectly proportional to the velocity of the Gibbs phase (elementary phase volume) change. 
This condition defines the general constant temperature systems. This condition leads 
to the canonical distribution as a stationary solution of the Liouville equations. For the 
linear friction, we derived the constant temperature systems. The general form of the 
non-potential forces is not considered in Ref. [12]. 
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• The canonical-dissipative systems are described in Ref. [9, 10]. These systems are the 
few-particle systems are defined by the non-potential forces = —dG(H)/dpi, where 
G(H) is a function of Hamiltonian H . The distribution functions are derived as solutions 
of the Fokker-Planck equation. Note that Fokker-Planck equation can be derived from 
the Liouville equation [11]. 

• The quantum few-particle systems with pure stationary states are suggested in Ref. [13, 
14]. The correspondent classical systems are not discussed. 

• The few-particle systems with the fractional phase space and non-Gaussian distributions 
are suggested in [15, 16]. Note that nondissipative systems with the usual phase space 
are dissipative systems in the fractional phase space [15, 16]. 

The analog of the first law of thermodynamics is connected with the variation of the mean 
value Hamiltonian 

U(x) = J if (q, p, x)p(q, p, x)d N qd N p (1) 

that has the form 

dU(x) = J 5 x Hpd N qd N p + J H5 x pd N qd N p. 

We can have the analog of the second laws of thermodynamics if the second term on the 
right-hand side can be represented in the form 

T(x)dS(x) = T(x)5 x J S N (p)pd N qd N p. 

In this case, we have the condition 

H5 x p = T(x)5 x (pS N (p)). 

This representation is realized if the distribution p can be represented as a function of Hamil- 
tonian p = p(H,x) such that we can write H = T(x)G(p). Obviously, we have the second 
requirement for this distribution: the distribution p satisfy the Liouville equation of the sys- 
tems. For these N-particle systems, we can use the analogs of the usual thermodynamics laws. 
Note that N is an arbitrary natural number since we do not use the condition > 1 or 
N — > oo. This allows us to use the suggested few-particle systems for the simulation schemes 
[17] for the molecular dynamics. 

In this paper we consider few-particle systems with distributions that are defined by Hamil- 
tonian and Liouville equation. We describe the few-particle systems that have some analog of 
the thermodynamic laws. These systems can be defined by the non-holonomic (non-integrable) 
constraint: the power of non-potential forces is directly proportional to the velocity of the el- 
ementary phase volume change. In the general case, the coefficient of this proportionality is 
determined by the Hamiltonian. The special constraint allows us to derive distributions for the 
system, even in far-from equilibrium states. The examples of these few-particle systems are the 
constant temperature systems [1, 2, 3, 4, 5, 6, 7, 8], the canonical-dissipative systems [9, 10], 
and the Fermi-Bose classical systems [9] . 
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In Sec. 2, we derive the analog of thermodynamic laws for the few-particle systems with 
the distributions that are defined by Hamiltonian. In Sec. 3, we consider the condition for 
the non-potential forces that allows us to use the analog of thermodynamic laws. We consider 
the wide class of few-particle systems with canonical Gibbs distribution and non-Gaussian 
distributions. In Sec. 4, we consider the non-holonomic constraint for few-particle systems. We 
formulate the proposition which allows us to derive the thermodynamic few-particle systems 
from the equations of few-particle system motion. The few-particle systems with the simple 
Hamiltonian and the simple non-potential forces are considered. Finally a short conclusion is 
given in Sec. 5. 



2 Thermodynamics Laws 
2.1 First Thermodynamics Law 

Let us consider the N-particle classical system in the Hamilton picture and denote the position 
of the ith particle by and its momentum by p,, where % — 1, N. The state of this system 
is described by the distribution function p = p(q, p, x, i). 

The mean value of the function /(q, p, x, t) is defined by the following equation 

f(x, t) = J /(q, p, x, t)p(q, p, x, t)d N qd N p. (2) 

Here, x = {x±, X2, ■■■x n \ are external parameters. The variation 5 X for this function can be 
defined by the relation 

<W(q, P, x,t) = 2^ £ dx k . (3) 

k=l k 

The first law of thermodynamics states that the internal energy U (x) may change because 
of (1) heat transfer SQ, and (2) work SA of thermodynamics forces Xk(x): 

n 

5A = J2Mx)dx k . (4) 



k=i 



The external parameters x here act as generalized coordinates. In the usual equilibrium ther- 
modynamics the work done does not entirely account for the change in the internal energy 
U(x). The internal energy also changes because of the transfer of heat, and so 

dU = 5Q- 5 A. (5) 

Since thermodynamic forces X k (x) are non-potential forces 

dX k (x) dXtix) 



dxi dx k 



(6) 



the amount of work 5 A depends on the path of transition from one state in parameters space 
to another. For this reason 5 A and 5Q, taken separately, are not total differentials. 
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Let us give statistical definitions of internal energy, thermodynamic forces and heat transfer 
for the few-particle systems in the mathematical expression of the analog of the first thermody- 
namics law. It would be natural to define the internal energy as the mean value of Hamiltonian 
H: 

U(x) = J if (q, p, x)p(q, p, x)d N qd N p. (7) 
It follows that the expression for the total differential has the form 

dU (x) = J S x H(q, p, x)p(q, p, x)d N qd N p + J if (q, p, x)S x p(q, p, x)d N qd N p. 

Using Eq. (3), we have 

dU (x) = J ^ ^fcp(q, P, x)d N qd N p + J if (q, p, x)S x p(q, p, x)d N qd N p. (8) 

In the first term on the right-hand side we can use the definition of phase density of the 
thermodynamic force 

Af(q,p,«) = 

ox k 

The thermodynamic force Xk(x) is the mean value of the phase density of the thermodynamics 
force 

X k (x) = J X d k (q, p, x)p(q, p, x)d N c L d N p. (9) 

Using this equation, we can prove relation (6). 

Analyzing these expressions, we see that the first term on the right-hand side of differential 
(8) answers for the work of thermodynamics forces (9), whereas the amount of the heat transfer 
is given by 

5Q = J if (q, p, x)5 xP (q, p, x)d N qd N p. (10) 

We see that the heat transfer term accounts for the change in the internal energy due not to 
the work of thermodynamics forces, but rather to the change in the distribution function cased 
by the external parameters x. 



2.2 Second Thermodynamics Law 

Now let us turn our attention to the analog of the second law for the few-particle systems. 
The second law of thermodynamics has the form 

SQ = T(x)dS(x). (11) 

This implies that there exists a function of state S(x) called entropy. The function T(x) acts 
as an integration factor. Let us prove that the law (11) follows from the statistical definition 
of SQ in Eq. (10). For Eq. (10), we take the distribution that is defined by the Hamiltonian, 
and show that Eq. (10) can be reduced to Eq. (11). 
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We can have the analog of the second laws of thermodynamics if the right-hand side of Eq. 
(10) can be represented in the form 

T(x)dS(x) = T(x)5 x J p(q,p,x)S N (p(ci,p,x))d N qd N p. 

This requirement can be written in an equivalent form 

H6 x p = T(x)S x (pS N (p)). 

This representation is realized if the distribution p can be represented as a function of Hamil- 
tonian 

p(q,p,x) = p(H(q,p,x),x). (12) 

Obviously, we have the second requirement for the distribution p, which must satisfies the 
Liouville equation. We assume that Eq. (12) can be solved in the form 

H = T(x)G(p), 

where G depends on the distribution p. The function T(x) is a function of the parameters 
x = {xi,x 2 , ...,x n }. As a result, we can rewrite Eq. (10) in the equivalent form 

5Q = j (T(x)G(p(q,p,x),x) +H (x))5 x p(q,p,x)d N qd N p. (13) 

The term with H (x), which is added into this equation, is equal to zero because of the nor- 
malization condition of distribution function p: 

H (x)S x J p(q,p,x)d N qd N p = H (x)S x l = 0. 

For canonical Gibbs distribution 

1 H(q,p,x) 

P = rvf eX P 



Z{x) * kT{x) ' 
where Z(x) is defined by 



u \ [ U{q,p,x N N 

(-)= j /exp-^---rf qd p, 



kT(x) 

we use G(p) and H Q (x) in the form 

H (x) = kT{x)lnZ{x) - kT(x), G(p) = -kln{Z{x)p). 
As a result, Eq. (13) can be rewritten in the form 

5Q = T(x)5 x J p(q, p, x)S N (p(q, p, x))d N qd N p, (14) 
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where the function Sn(p) is defined by 

d(pS N ( P )) 
dp 



G(p) + H (x)/T(x). (15) 



We see that the expression for SQ is integrable. If we take 1/T(x) for the integration factor, 
thus identifying T(x) with the analog of absolute temperature, then, using Eqs. (11) and (14), 
we can give the statistical definition of entropy: 

S(x) = J p(q, p, x)S N (p(q, p, x))d N qd N p + S . (16) 

Here So is the contribution to the entropy which does not depend on the variables x, but may 
depend on the number of particles N in the system. As a result, the expression for entropy is 
equivalent to the mean value of phase density function S d (q,p,x) = S(p(q,p,x)) + Sq. Here 
S d is a function of dynamic variables q, p, and the parameters x = {xi, X2, ...,x n }. The number 
N is an arbitrary natural number since we do not use the condition iV ^> 1 or N — > oo. Note 
that in the usual equilibrium thermodynamics the function T(x) is a mean value of kinetic 
energy. In the suggested thermodynamics of few-particle systems T(x) is the usual function of 
the external parameters x = {xi,x 2 , ...,x n }. 



2.3 Thermodynamic Few-Particle Systems 

Let us define the special class of the few-particle systems with distribution functions that are 
completely characterized by the Hamiltonian. These distributions must satisfy the Liouville 
equation for the few-particle system. 
Definition A few-particle system 

<k\j_ = dH_ dpi_ = _dH_ p ( n ) 
dt dpi dt dqi 1 

will be called a thermodynamic few-particle system if the following conditions are satisfied: 

(1) the distribution function p is determined by the Hamiltonian, i.e., p(q,p,x) can be written 
in the form 

P(q, P, x) = p(H(q, p, x),x), (17) 
where x is a set of external parameters, and 

p(q,P,z)>0, J p(q,p,x)d N qd N p = 1; 

(2) the distribution function p satisfies the Liouville equation 

i+|-(^)+|-(<-g+^M 

(3) the number of particles N is an arbitrary natural number. 
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Here and later we mean the sum on the repeated index i from 1 to N. 
Examples of the thermodynamic few-particle systems: 

(1) The constant temperature systems [1, 2, 3, 4, 5, 6, 7, 8] that have the canonical distribution. 
In general, these systems can be defined by the non-holonomic constraint, which is suggested 
in Ref. [12]. 

(2) The classical system with the Breit-Wigner distribution function that is defined by 

^"-'^ M^FiMVf ' (18) 

(3) The classical Fermi-Bose canonical-dissipative systems [9] that are defined by the distribu- 
tion functions in the form 

3 Distribution for Thermodynamic Few- Particle Systems 
3.1 Formulation of the Results 

Let us consider the few-particle systems which are defined by the equations 

+ FS n) , (20) 



dcii dH dpi dH i ^ (n) 



dt dp/ dt dq, 
where % = 1, N. The power of non-potential forces is defined by 

a tt 

V( q ,p,x) = F^^±. (21) 

If the power of the non-potential forces is equal to zero (V = 0) and dH/dt = 0, then few- 
particle system is called a conservative system. The velocity of an elementary phase volume 
change fl is defined by the equation 

_ dFj d 2 H _ dF^ 
q,p,x dpi dqidpi dpi 

We use the following notations for the scalar product 

dAj _ ^, dA xi dA yi dAzj \ 
d&i ~{^da xi da yi da zi '' 

The aim of this section is to prove the following result. 
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Proposition 1. If the non-potential forces of the few-particle system (20) satisfy the 
constraint condition 

3F, W _ ( „>£W 



a P , -^ H ' x)F 'Wr ' (23) 

then this system is a thermodynamic few-particle system with the distribution function 

p(q, p, x) = exp —B(H(q, p, x),x), (24) 

Zj\X) 

where the function B(H,x) is defined by the equation dB(H,x)/dH = /3(H,x), and 

Z(x) = J exp[-B(H(q,p,x),x)]d N qd N p 

Obviously, we consider the distribution functions (24) and the function B = B(H,x) such 
that p > 0, and Z(x) < oo. Note that condition (23) means that the velocity of the elementary 
phase volume change f2 is directly proportional to the power V of non-potential forces 
of the few-particle system (20) and coefficient of this proportionality is a function (3(H,x) of 
Hamiltonian H, i.e., 

fi(q, p, x) - (3(H, x)V(q, p, x) = 0. (25) 

Note that any few-particle system with the non-holonomic constraint (25) or (23) is a 
thermodynamic few-particle system. Solving the Liouville equation with the non-holonomic 
constraint (23), we can obtain the distributions that are defined by the Hamiltonian. 

3.2 Proof of Proposition 1 

Let us consider the Liouville equation for the few-particle distribution function p = p(q, p, x, t). 
This distribution function p(q, p, x, t) express a probability that a phase space point (q, p) will 
appear. The Liouville equation for this few-particle system 



t + l-M+^M^O (26) 



expresses the conservation of probability in the phase space. Here we use 

dH „ OH , n) 

K = F- = h F ( ' 

dpi dqt 

Using a total time derivative along the phase space trajectory by 

d d d d 

dt dt 1 dqi % ' dp^ 

we can rewrite Eq. (26) in the form: 
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where the omega function is defined by Eq. (22). In classical mechanics of Hamiltonian systems 
the right-hand side of the Liouville equation (28) is zero, and the distribution function does not 
change with time. For the N-particle systems (20), the omega function (22) does not vanish. 
For this system, the omega function is defined by (22). For the thermodynamic few-particle 
systems, this function is defined by the constraint (23) in the form 

n = ^x)Fi"^. (29) 

In this case, the Liouville equation has the form 

^Il = . mx)¥ ^ p . (30) 

Let us consider the total time derivative of the Hamiltonian. Using equations of motion (20), 
we have 



dB__dE_ dE_dE_ / dH F (n)\ dH _ OH {n )dH_ 



If dH/dt = 0, then the power V of non-potential forces is equal to the total time derivative of 
the Hamiltonian 

1 dpi dt 
Therefore Eq. (30) can be written in the form 

dt dt 
If f3(H,x) is an integrable function, then this function can be represented as a derivative 

,33) 

In this case, we can write Eq. (32) in the form 

din p(q, p, x) dB(H,x) 
di ~~ dt ' 

As a result, we have the following solution of the Liouville equation 



(34) 



p(q, p, x) = exp -B(H(q, p, x),x). (35) 

The function Z(x) is defined by the normalization condition. It is easy to see that the dis- 
tribution function of the N-particle system is determined by the Hamiltonian. Therefore, this 
system is a thermodynamic few-particle system. Note that iV is an arbitrary natural number 
since we do not use the condition iV 3> 1 or N — > oo. 
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3.3 Few-Particle Systems with Canonical Distributions 

In this section, we consider the thermodynamic few-particle system that is described by canon- 
ical distribution [12]. These few-particle systems are defined by the simple function f3(H,x) = 
3N/3(x) in the non-holonomic constraint (25). 

Corollary 1. // velocity of the elementary phase volume change fl is directly proportional 
to the power of non-potential forces V , then we have the usual canonical Gibbs distribution as 
a solution of the Liouville equation. 

In other words, the few-particle system with the non-holonomic constraint Q = (3(x)V can 
have the canonical Gibbs distribution 



p(q, p, x) = exp p(x) (F(x) - H(q, p, 



x 



as a solution of the Liouville equation. Here the coefficient (3(x) does not depend on (q, p, t), 
i.e., df3(x)/dt = 0. Proof of this corollary is considered in Ref. [12]. 
Using Eq. (21), we get the Liouville equation in the form 

MS^) = _ /?WF (»»^ p . (36 ) 

The total time derivative for the Hamiltonian is defined by Eq. (31) in the form 

dH = dH_ + F (n)dH_ 
dt dt 1 dpi 

If dH /dt = 0, then the energy change is equal to the power V of the non-potential forces F\ n \ 
Therefore the Liouville equation can be rewritten in the form 

dlnp(q,p,x) Q( ^ dH(q,p,x) _ 
dt +IJ[X) dt " U - 

Since coefficient f3(x) is a constant (d(3(x)/dt = 0), we have 

j t (lnp(q., p, x) + P(x)H(q, p, xfj = 0, 

i.e., the value (Inp + f3H) is a constant along of the trajectory of the system in 6N-dimensional 
phase space. Let us denote this constant value by (3(x)J 7 (x). Then we have 

lnp(q, p, x) + P(x)H(q, p, x) = / 9(x)J r (a;), 

where dT{x)jdt = 0. As a result, we get a canonical distribution function 

p(q, p, x) = exp p(x) (F{x) - H(q, p, xj) . 

The value J-{x) is defined by the normalization condition. Therefore the distribution of this 
few-particle system is a canonical distribution. 
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3.4 Non-Canonical Distributions for Few-Particle Systems 

In Sec 3.3, we consider (3(H,x) = f3(x). Let us consider the linear function (3 = (3(H,x). 

Corollary 2. The linear function /3(H,x) in the form 

(3(H,x) = (3 1 (x)+(3 2 (x)H 
leads to the following non-canonical distribution function 

p(q,p,x) = -^^exp-^P^H +^p 2 (x)H 2 y 

The proof of this proposition can be directly derived from Eqs. (35) and (33). 

The well known non-Gaussian distribution is the Breit-Wigner distribution. This distri- 
bution has a probability density function in the form p(x) = l/7r(l + x 2 ). The Breit-Wigner 
distribution is also known in statistics as Cauchy distribution. The Breit-Wigner distribution 
is a generalized form originally introduced [18] to describe the cross-section of resonant nuclear 
scattering in the form 

" Wq ' p ' x)) = (g(q, P ,s)-g)» + (r/2)' - (37) 

This distribution can be derived from the transition probability of a resonant state with known 
lifetime [19, 20, 21]. The second non-Gaussian distribution, which is considered in this sec- 
tion, is classical Fermi-Bose distribution that was suggested by Ebeling in Refs. [9, 10]. This 
distribution has the form 

Corollary 3. If the function (3(H,x) of the non-holonomic constrain is defined by 

^= ( H -Ei + %,2r (39) 

then we have thermodynamic few-particle systems with the Breit-Wigner distribution (37). 
Corollary 4. If the function j3(H,x) of the non-holonomic constrain has the form 

P(H,X) = — m R( ,rr , (40) 

1 + a exp (3(x)H 

then we have thermodynamic few-particle systems with classical Fermi-Bose distribution (38). 

Note that Ebeling derives the Fermi-Bose distribution function as a solution of the Fokker- 
Planck equation. It is known that Fokker-Planck equation can be derived from the Liouville 
equation [11]. We derive classical Fermi-Bose distribution as a solution of the Liouville equation. 
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If the non-potential forces are determined by the Hamiltonian 

F< B > = -dOW/dpi, (41) 

then we have the thermodynamic few-particle systems, which are considered in Refs. [9, 10]. 
These systems are called canonical dissipative systems. 
Let us assume that Eq. (17) can be solved in the form 

H = T(x)G(p), (42) 

where G depends on the distribution p. The function T(x) is a function of the parameters x. 
In this case, the function f3(H,x) is a composite function 

C(p) = -(3(T(x)G(p),x). (43) 

This function can be defined by 

As a result, we have the Liouville equation for the few-particle system in the form 

g = CMP. (45) 

This equation is a nonlinear equation. For example, the classical Fermi-Bose systems [9] have 
the function in the form 

C(p) = -P(x)(p-sp 2 ). (46) 

Note that the nonlinear evolution of statistical systems is considered in Refs. [22, 23, 24, 25, 
26, 27, 28, 29]. 



4 Non-Holonomic Constraint 



4.1 Formulation of the Result 



In this section, we formulate the proposition, which allows us to derive the thermodynamic 
few-particle systems from any equations of motion of N-particle systems. The aim of this sec- 
tion is to prove the following result. 



Proposition 2. For any few-particle system, which is defined by the equation 

+ F< B \ i = l,..,JV, (47) 



dq { OH dpi OH : Wn) 



dt d~pi dt <9q. 
there exists a thermodynamic few-particle system that is defined by the equations 

d(\i dH dp 



dt dpi ' dt 



= Fr w , (48) 
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and the distribution (24), where the non-potential forces Ff ew are defined by 

F new _ PfcPfc^j ~ PjPj ( dH (ra) N P^Qj dH 
1 " P fc P fc V dqj ^ i ) PfcPfc dpj ■ ^ > 

The vectors Pj and Qj are defined by the equations 

= dmx) dHdH F(n) + 9Ff9ff + 

<7// <9pi c7p ; J ' <9p, <9pj 

+ftH,x)F™ (50) 

J dpidpj opidpj 

Q = dp(H,x) dHdH F(n) | dF^dH^ 
dH dqi dpj j ' dqi dpj 

+P (H, x )Ff/^--^. (51) 
dqidpj dqidpj 



and 



Here we use the following notations 

N 

ajhjCj = a, ^^{b X jC X j + b y jC y j + b Z jC Z j). 

3=1 

Note that the forces that are defined by Eqs. (49), (50) and (51) satisfy the non-holonomic 
constraint (25), i.e., 

d¥f w | d 2 H [3<H,x)F newdH (52) 
dpj dqjdpj ' j dpj 

where we use the omega function in the form (22). 
4.2 Proof of Proposition 2. 

Let us prove Eq. (49). Let us consider the N-particle classical system in the Hamilton picture. 
Denote the position of the ith particle by and its momentum by p^. Suppose that the system 
is subjected to a non-holonomic (non-integrable) constraint in the form 

/(q,p,x)=0. (53) 

Differentiation of Eq. (53) with respect to time gives a relation 

P,(q, p,x)^ + Q,(q, p,x)^ = 0, (54) 

where the functions Pj and Qj are defined by the equations 

P;(q, p, x) = J^, Q 4 (q, p, x) = J^. (55) 
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An unconstrained motion of the ith particle, where i — 1, N, is described by the equations 

^ = K>, ^ = F 4 , (56) 
dt dt K 1 

where Fj is a resulting force, which acts on the ith particle. 

The unconstrained motion gives a trajectory which leaves the constraint hypersurface (53). 
The constraint forces Rj must be added to the equation of motion to prevent the deviation 
from the constraint hypersurface: 

^ = K l? ^ = F, + R, (57) 
dt ' dt K ' 

The constraint force Rj for the non-holonomic constraint is proportional to the [30] : 

R = APi, (58) 

where the coefficient A of the constraint force term is an undetermined Lagrangian multiplier. 
For the non-holonomic constraint (53), the equations of motion (56) are modified as 

^ = Ki , = F, + AP,. (59) 

dt dt K J 

The Lagrangian coefficient A is determined by Eq. (54). Substituting Eq. (57) into Eq. 
(54), we get 

Pi(Fi + APi) + QjKj = 0. (60) 
Therefore the Lagrange multiplier A is equal to 

P t F t + Q t K t 

A = B-5 • ( 61 ) 

As a result, we obtain the following equations 

dt dt P fc P fc v ; 

These equations we can rewrite in the form (56) 

^ = K l5 ^ = F»» (63) 
at at 



with the new forces 



PfcPfe^jj P«Pj T-\ PjQj 



op F , - (64) 

In general, the forces F^ ew are non-potentials forces (see examples in Ref. [12]). 

Equations (62) are equations of the holonomic system. For any trajectory of the system 
in the phase space, we have / = const. If initial values qfc(0) and Pfc(0) satisfy the constraint 
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condition /(q(0), p(0), x) = 0, then solution of Eqs. (62) and (64) is a motion of the non- 
holonomic system. 

Let us prove Eqs. (50) and (51). In order to prove these equations, we consider the few- 
particle system (56) with 

K, = jj» F, = -^ + F<">, (65) 
dpi <9q; 

and the special form of the non-holonomic constraint (53). Let us assume the following con- 
straint: the velocity of the elementary phase volume change f2(q, p, x) is directly proportional 
to the power V(<i, p, x) of the non-potential forces, i.e., 

Q(q,p,x)=P(H,x)V(q,p,x), (66) 

where /3(H,x) depends on the Hamiltonian H. Therefore the system is subjected to non- 
holonomic (non-integrable) constraint (53) in the form 

/(q, p, x) = (3(H, x)V(q, p, x) - Q(q, p, x) = 0. (67) 

This constraint is a generalization of the condition which is suggested in Ref. [12]. The power 
V of the non-potential forces F- n ^ is defined by Eq. (21). The function Vt is defined by Eq. 
(22). As a result, we have Eq. (67) for the non-potential forces in the form 



[3(H,x)F*f ) ^-^-^ = 0. 



dH <9F$ n) 
d Pj dpj 



The functions Pj and Qj for this constraint can be find by differentiation of constraint. Differ- 
entiation of the function /(q, p, x) with respect to Pi gives 



Op, dpi v - 3 dpj > dpi dpj 



This expression leads us to Eq. (50). Differentiation of the function /(q, p, t) with respect to 
q* gives 



w™4 = £ = £(M*rt ,! B 

This expression leads to Eq. (51) 



dqi dqS^^ H,X ^ Fj dp^ dc[i dpj 



4.3 Few-Particle Systems with Minimal Constraint 

Let us consider the simple constraints for to realize the classical systems with canonical and non- 
canonical distributions. Let us consider few-particle system, which is defined by the simplest 
form of the Hamiltonian 

H( (l ,p,x) = ^ + U(q,x), (68) 
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and the non-potential forces 

F\ n) = (69) 

Here p 2 = ££i p 2 , and N is a number of particles. For the minimal constraint models, the 
non-holonomic constraint is defined by the equation 

f(q,p,x)=P(H,x)^--3N = 0. (70) 

777* 

The phase space gradients (50) and (51) of this constraint are represented in the form 

_ ( df3(H,x) p 2 \ 2p, d/3(H,x) dH 

OH 2m +IJ[H,X) ) m> ^ ~ OH dq- 

The non-potential forces of the minimal constraint models have the form 

F ne W = P% ~ PiPj dU , PjVj 9 P( H > X ) dU 

1 p 2 9 qj 2p 2 ([p 2 /2m)d(3(H,x)/dH] +0(H,x)) dH 9q/ 

It is easy to see that all minimal constraint models have the potential forces. 
Note that the minimal Gaussian constraint model is characterized by 

d(3(H,x) 

In this case, we have the non-potential forces in the form 



p 2 dq, ■ 

The few-particle systems are the constant temperature systems that are considered in Refs. 
[1, 2, 3, 4, 5, 6, 7, 8, 12]. 

4.4 Few-Particle Systems with Minimal Gaussian Constraint 

Let us consider the N-particle system with the Hamiltonian (68), the function /3(H,x) = 
3N/kT(x), and the linear friction force (69). Substituting Eq. (69) into Eqs. (21) and (22), we 
get the power V and the omega function Q: 

V = -— p 2 , n = -3 7 iV. 
m 

In this case, the non-holonomic constraint has the form 

= kT(x), (71) 

m 

i.e., the kinetic energy of the system must be a constant. Note that Eq. (71) has not the 
friction parameter 7. 
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For the few-particle system with friction (69) and non-holonomic constraint (71), we have 
the following equations of motion 

dqi pi dpi dU df 
at m at oqi dpi 

where the function / = /(q, p) is defined by 

f(^p) = \{p 2 -mkT(x)) : /(q,p)=0. (73) 

Equation (72) and condition (73) define 6N+1 variables (q, p, A). Let us find the Lagrange 
multiplier A. Substituting Eq. (73) into Eq. (72), we get 

^ = -9q" + (A - 7)P - (74) 



Using df/dt — in the form 



P*i = (75) 



and substituting Eq. (74) into Eq. (75), we get the Lagrange multiplier A in the form 

_1 dU_ 

A ~ mkT{x) Pj dq~ 4 

As a result, we have the holonomic system that is defined by the equations 



rfqi _ pi dpi_ _ 1 dU _ dU_ . . 

dt m' dt mkT{x) 13 d<\j dc^ 

For the few-particle system (76), condition (71) is satisfied. If the time evolution of the few- 
particle system is defined by Eq. (76), then we have the canonical distribution function in the 
form 

P(q ' P '" ) = ^ eXp -^^- (77) 

where Z(x) is defined by the normalization condition. For example, the few-particle system 
with the forces 

2 ( \ 

F, = -mu 2 (x)q t + ^g-p^p^.) (78) 
has canonical distribution (77) of the linear harmonic oscillator with Hamiltonian 

tt, \ P 2 ™ 2 (i)q 2 
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5 Conclusion 



In this paper we derive the extension of the statistical thermodynamics to the wide class of 
few-particle systems. We consider few-particle systems with distributions that are defined 
by Hamiltonian and Liouville equation. These systems are described by the non-holonomic 
(non-integrable) constraint [12]: the velocity of the elementary phase volume change is directly 
proportional to the power of non-potential forces. In the general case, the coefficient of this 
proportionality is defined by Hamiltonian. This constraint allows us to derive the distribution 
function of the few-particle system, even in far-from equilibrium states. The few-particle sys- 
tems that have some analog of the thermodynamic laws is characterized by the distribution 
functions that are determined by the Hamiltonian. The examples of these few-particle systems 
are the constant temperature systems [1, 2, 3, 4, 5, 6, 7, 8], the canonical-dissipative systems 
[9, 10], and the Fermi-Bose classical systems [9]. For the few-particle systems, we can use the 
analogs of the usual thermodynamics laws. Note that the number of particles is an arbitrary 
natural number since we do not use the condition N 1 or N — > oo. This allows one to use 
the suggested few-particle systems for the simulation [17] for the molecular dynamics. 

The quantization of the evolution equations for non-Hamiltonian and dissipative systems 
was suggested in Refs. [31, 32]. Using this quantization it is easy to derive the quantum analog 
of few-particle systems that leads to some analog of the thermodynamic laws. We can derive the 
canonical and non-canonical statistical operators [33] that are determined by the Hamiltonian 
[13, 14]. The suggested few-particle systems can be generalized by the quantization method 
that is considered in Refs. [31, 32]. 
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